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DEFINITION OF SYMBOLS 


Symbol 


k 


1 


m, m^ 
2r 
t 


Q i 


y., z. 
J V i 


C 


E 

E 

n 

F 

I 

P 

Q i 

S 2 

T 

T 

n 

U 


Definition 

Spring constant of a torsional spring, inch -pounds . 
Spring constant of a linear extensional spring, pounds. 
Length of link, inches. 

Mass, pounds. 

Length of major axis of an ellipse. 

Time, seconds. 

Generalized coordinates of model. 

Rotated set of coordinates . 

Constant of Euler load = . 

Total energy of the system, inch-pounds. 

Total energy of the nondimens ional system. 

p 

Magnitude of the nondimensional load, — — 

cr 


Magnitude of the nondimensional impulse, At . 
Magnitude of load, pounds. 

Generalized forces, inch-pounds. 

, , , ml 2 co 2 7 -n/ 41 

Dimensionless constant, : = — . 

k 4 

Kinetic energy, inch-pounds. 

Nondimensional kinetic energy. 

Potential energy of internal forces or strain energy, inch- 
pounds . 
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Symbol 
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6, A 

e. 
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<P i 
<P 2 
4 > 

CO 

S2 

S2 

n 


Definition 

Nondimensional potential energy of internal forces. 

Total potential energy, inch-pounds. 

Nondimensional total potential energy. 

Arbitrary constant that softens the spring. 

Incremental change in variables . 

Imperfection or eccentricity of links, radians. 

Changed coordinates where cp. = 

1 ^/2 

Nondimensional time, cot. 

Angular displacement of lower link from the vertical, radians. 
Angular displacement of upper link from the vertical, radians. 
Angular rotation of y., z. axes, degrees. 

Angular frequency, radians/second. 

Potential energy of external load, inch-pounds’. 

Nondimensional potential energy of external load. 


Superscripts Definition 

Derivative of displacement with respect to t. 

Second derivative of displacement with respect to t. 
! Derivative of displacement with respect to r. 
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DEFINITION OF SYMBOLS (Concluded) 


Superscripts Definition 

" Second derivative of displacement with respect to t. 

a Initial velocity after impulse . 

~ Evaluated in the reference state. 

0 Degree . 

Subscripts Definition 

b Initial energy, 

cr Critical value . 

D Dynamic load, 

f Final energy. 

i, j, k, 1 Index notation with range of 1 to 2 except as noted, 

n Nondimens ional value. 

0 Initial condition. 

S Static load. 
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TECHNICAL MEMORANDUM X-64531 


STABILITY OF A TWO-DEGREE-OF- FREEDOM SYSTEM 
WITH STATIC AND DYNAMIC LOADING 

SECTION I. SUMMARY 


Energy methods are applied to determine the stability of a nondimen- 
sional two-degree -of -freedom column model under static and impulsive 
loading. A soft spring is used to create instability, and provision is made for 
initial imperfections. No dissipation force is considered. The principle of 
minimum potential energy is used to determine stability of the linear system 
at the critical static buckling load. A two-dimensional potential energy 
contour plot locates the saddles of the nonlinear deflected system and the 
critical potential energy of the saddle is determined. Stability at the saddle 
without a static load is investigated. 

The initial kinetic energy to reach the saddles with and without static 
loading is provided by a pure impulse ( initial velocities of the system) or a 
square-wave impulsive load with an initial imperfection. Solution curves are 
determined for the various initial conditions. A numerical integration is used 
to determine the system response and the excess initial kinetic energy required 
to cause collapse is determined for various conditions . 

A relationship is shown between the pure impulse and the impulsive 
load with imperfections. The difference in stability criterion at the initial 
static critical buckling load and at the saddle of the deflected system without a 
static load is noted. The effects of the various system parameters are dis- 
cussed in the conclusions. 

This is an initial investigation of a simplified system. Investigation of 
more complex systems may lead to application of the methods to physical 
systems . 



SECTION!!. INTRODUCTION 


A. Buckling 

In a broad sense, the theory of buckling is the theory of the stability 
of equilibrium of mechanical systems . As the load on a mechanical system 
increases, the equilibrium states trace paths in a configuration space. At 
some load a "bifurcation point" or fork may be reached, and different paths 
may be traced from this point. Exceeding this critical (or Euler) load 
results in buckling, and the initial postbuckling behavior of the system depends 
on the stability or instability of equilibrium at the critical buckling load. 

Most of the phenomena associated with instability and buckling can be 
observed in simple systems. The one -degree -of -freedom system depicted in 
Figure la illustrates this phenomenon.* It is composed of two rigid bars joined 
by a frictionless pin, and constrained by a linear extensional spring (k t ) and 
a torsional spring (k) . The top is constrained to move vertically while the 
bottom is pinned to a fixed support by a frictionless pin. 

B. Bifurcation of Equilibrium States 

Figure lb depicts the deflection of the system as a function of a 
nondimensional load. At the bifurcation point (P = P ) the equilibrium paths 

branch into the curves 1A or IB. The horizontal branch depicts the classical 
theory. 


If the extensional spring (ki) is removed, the equilibrium path is 
1A, which represents a stable configuration. This is typical of the actual 
postbuckling behavior of large bars. If the extensional spring is retained and 
the torsional spring (k) is removed, the equilibrium path is IB, which 
represents unstable configurations. In some types of structures, particularly 
shells, an instantaneous displacement of considerable magnitude can occur 
as the structure moves from one point on the equilibrium path to another point 
on the path. This is known as "snap -through" and is shown in Figure 2. 

Imperfection-sensitivity and its effect on static buckling has been 
studied intensively in recent years. The curves of Figure 3 are based on 
the studies of Koiter [i ] , where the solid lines represent the states of 


* This model was developed by Dr. G. A. Wempner, University of Alabama 
in Huntsville . 
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a b 

Figure i. Buckling of a column. 



equilibrium of perfect systems, and 
the dotted lines represent the states of 
equilibrium of imperfect systems. 

Figure 3a represents stable equilibrium 
as the load must increase with deflection. 
Figure 3b represents unstable equilib- 
rium as the load decreases with deflec- 
tion. Figure 3b has greater imperfec- 
tion sensitivity. 

C. Purpose of the Study 


Figure 2. Snap-through. The effects of static and dynamic 

loading on a simplified nondimensional 

two-degree -of -freedom system are examined to determine whether the methods 
are applicable to stability analyses. The study contains four parts: 


1. A linear analysis to determine the critical conditions for initial 
buckling. 

2. A nonlinear analysis of the potential energy to ascertain the stability 
of the equilibrium states. 
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DISPLACEMENT 
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Figure 3. Postbuckling behavior and imperfections. 

3. An analysis of the dynamic loading that will cause unrestricted 
deflections . 


4. An analysis of the effect of initial imperfections. 

The principle of minimum potential energy and a Taylor series 
expansion about the equilibrium state are used in the stability analysis . No 
dissipation force is considered. Lagrange's equations of motion for 
conservative systems are used in the dynamic analysis with an applied static 
and/or dynamic load. The initial kinetic energy that will equal the critical 
potential of the saddle is determined. A numerical integration is then used 
to determine the excess initial kinetic energy that will result in collapse of 
the system. 


SECTION III. THE MODEL 
A. Description of the Model 

The two-degree -of -freedom model used in the analyses is shown in 
Figure 4. it is a double pendulum composed of two rigid weightless bars of 
equal length, 1, which carry two masses such that m t = 2m and m 2 = m. 

The pins on which the bars rotate are frictionless. The generalized coordinates 
are angles (p t and <p 2 measured from the vertical. An imperfection, a 
small angle ( e ) , is provided for each link to determine the effect of imper- 
fections under impulse loading. The restoring force consists of a nonlinear 
torsional spring at each pin such that the restoring force is of the form 
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k(q>- /3(p 3 ) 

with /? an arbitrary constant. The soft spring is to produce conditions 

similar to snap-through. A force is applied at the free end of the linkage; 

this is a conservative static force (P ) and/or a dynamic force ( P ) where 

b D 

P is a function of time. The force of gravity on the masses is neglected. 

The kinetic energy, T, the external potential energy, J2, the 
internal potential energy, U, and the generalized forces, Q* and Q 2 , of 
the system shown in Figure 4 are 


T 

(2 

U 



ml 2 [3<pi + 2<p 1 <p 2 cos( cp 2 + e 2 - <p x - e t ) + <j? 2 2 ] 
1 [cos^ + £l ) + eos(<p 2 + e 2 ) - 2] , 

J2L _ + LS£l Z £l l l _ a iV'L.ZJPl)! 

2 P 4 2 p 4 


9 


9 


(la) 

(lb) 

(lc) 
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Qi = P D 1 sin ( (Pi + € i) 

(id) 

Q 2 = P-pl sin( cp 2 + e 2 ) 

(ie) 


The total potential energy of the system is 

V = U + £2 (2) 

B. Requirements for Stable Equilibrium 

The principle of minimum potential energy states that a conservative 
holonomic system is in a configuration of stable equilibrium if, and only if, the 
potential energy is a relative minimum [2] . To examine the requirements for 
stable equilibrium, let q. denote the generalized coordinates of the system, 

with i = 1, 2. Then the change in potential energy is 
AV = V(q.+ 6q.) - V(q.) . 

Expansion of this in a Taylor series can be written as 


AV = SV + 6 2 V + d^V + <5 4 V + R (3) 

where oV is known as the first variation, S 2 V as the second variation, 
etc. 


The system has a stationary value and is in equilibrium if 6V= 0. 
This stationary value can be a maximum, a minimum, or a saddle, and the 
state of stability is determined by the higher order terms of equation (3) . 
Since the potential energy is a relative minimum for stable equilibrium, 
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AV > 0 is necessary and sufficient for stability [i] . A further necessary 
condition for stability [i] is 

S 2 V = 0 . (4) 


C. Equations of Equilibrium 

For the two-degree -of -freedom system shown in Figure 4, the first 
variation is 



<Vi + 


8V 

dcp 2 


d(p 2 


Since 6<p i and 6cp 2 are independent, the equilibrium conditions can be 
expressed as 


9 V 


= 0 


5 


and the equations of equilibrium are 


-Pgl sin(^! + Ci) + k{2^ 1 - <p 2 - [ <p ± 3 - (cp 2 ~ <P \) 3 ]} = 0 , (5a) 

-P g l sin(<p 2 + e 2 ) +k[^ 2 “ <Pi ~ 0(<?2 " <Pi) 3 ] =0 • (5b) 


D. Critical Buckling Load 


The second variation, 6 2 V, is a positive definite quadratic when 

P < P . With 6V = 0, the second variation is 
cr 


6 2 V = 


1 / 9 2 V 

2 *. \ 


*P i 


+ 2 


9 2 V 

9(^1 9 <P2 


5( Pl 6( P2 + 


9 2 V 

~<P2 



The right side can be expressed in matrix form as X $ , or 
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5 2 V= [6cpi 6<p 2 ] 


a 2 v a 2 v 


6«pi 

9<Pl Z &<Pi 9<^2 


9 2 V 9 2 V 


6<p 2 

d(pid(p 2 9 <pi 



( 6 ) 


For other than trivial conditions, the character of equation (6) is determined 
by the matrix X . At the critical buckling load, the matrix X becomes 
positive semidefinite [2] and 

det X = 0 . 

In determining the critical buckling load of the perfect system, e. is 

neglected. The higher order terms of <p. are also neglected since near 
cpi = cp 2 = 0 these are much smaller than < ■p i and cp 2 . Substitution of 
equations (lb), (ie), and (2) into equation ( 6) and setting det X = 0 yields 


-PI + 2k -k 

-k -P 1 + k 


= 0 . 


Solution of this gives 


P 


/ 3 ±n/T 

V 2 


) 


k 

1 




and the lower value is the critical or Euler buckling load 


P 

cr 


/ 3 - nTT 
\ 2 


) 


k 

1 


(7) 
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E. Initial Buckling Mode 

To determine the initial buckling mode, P is substituted for P 

cr 

in equation (5a) or (5b) , again neglecting higher orders of <p.. Equation 
( 5a) then gives the initial buckling mode as 



F. Equations of Motion 

Since the system of Figure 4 is holonomic, the nonlinear equations of 
motion are given by Lagrange's equations [3] as 

_d_ / 8T j _ 3T + _8V =Q 

dt \'*i m d(f> i d(f>i i 


Substitution of equation ( i) into this gives 

ml 2 [3pi+ <p2 cos ( cp2 + e 2 “ <Pi “ e i) " cp2 2 sin(<p 2 + ^2 ~ <Pi ~ ei)l 

(9a) 

+ k (2 <p x - <p 2 - @[<pi - (<p 2 - <?i) 3 ]} = (P g + P D ) 1 sin (<Pi + e*) , 

ml 2 [‘<p 2 + t Pt cos ( (p2 + £2 “ <Pi ~ ^i) + (f>t sin(<p 2 + ~ <Pt " <4)] 

(9b) 

+ k[(p 2 - <p i - P ( <P2 ~ <P \) 3 ] = ( p g + p D ) 1 sin( (p 2 + €2) 

To determine the linear frequency equation with P^ = 0, assume 
solutions of equation ( 9) in the form of 

c p = A. sin(o;t + A ) 

with (pi « 1 and <p 2 « 1, which yields 

( -3ml 2 w 2 + 2k - P„l) Ai + ( -ml 2 w 2 - k) A 2 = 0 , 

O 

( -ml 2 o> 2 - k) A* + ( -ml 2 w 2 + k - Pgl) A 2 = 0 . 
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This has a nontrivial solution only if the determinant of the coefficients of 
and A 2 vanish. Let 


ml 2 gj 2 

k 


= S 2 



= N . 


The frequency equation is then 


2S 4 + (4N - 7) S 2 + (N 2 - 3N + 1) = 0 


( 10 ) 


If P = 0, this reduces to 
b 

2S 4 - 7S 2 + 1=0 


and 

S 2 = 7±n/~41 


The lowest natural frequency of the linear system is then 


S 2 


7 - n/~41 
4 


( 11 ) 


To determine the critical load of the linear system, let co - 0. Then 
S 2 = 0 and equation ( 16) becomes 

N 2 - 3N + 1 = 0 , 


and 
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N = 


P 1 

S 3 ±\l~5 


The critical or buckling load is the lower value 


cr 


/ 3 - n/~ 5 \ _k 

V 2 / 1 


which agrees with equation ( 7) . 

G. Nondimensional Equations 

To eliminate certain physical parameters, the equations of motion 
are changed to a nondimensional form by using the dimensionless quantities: 


S 2 = 


ml 2 co 2 


(12a) 


F S P 

cr 


(12b) 


T = Wt 


(12c) 


With cp' = wcp and <p" = c o 2 lp , the nondimensional equations of motion are 


S 2 [3cp i " + (p 2 " cos(cp 2 + e 2 - cpt - €i) - (<p 2 ') 2 sin(<p 2 + e 2 - <pi - e^] 


+ { 2 cpi - <p 2 ~ P[<p? - ((Pi ~ <Pi) 3 ]} = CF sin( (pi + ej) , 


(13a) 


S 2 |> 2 " + cp t " cos((p 2 + e 2 - (pi - ej) + (cp i ') 2 sin(^ 2 + *2 ~ <Pi~ ^l)] 


+ {<P2 - <Pi - <p2 - <Pi) 3 ) = CF sin( <p 2 + e 2 ) , 


(13b) 


where S 2 = — — C = and F is defined as F^ and/or F„. 

4 2 D S 
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The nondimensional equations of equilibrium are 


2 <p x - cp 2 ~P[(pt ~ (<P 2 ~ <Pi) S ] ~ CF g s in(<p i + ei) = 0 , (14a) 

CP2 - (Pi - (H<P2 - <Pi) 3 " CF g sin ( <P2 + e 2> = 0 • ( 14b ) 

The nondimensional energy equation is 
S 2 

E n = — [3(^i») 2 + 2 ^ 1 * <p 2 ' cos(<p 2 + € 2 - (pi - ei) + (<p 2 f ) 2 ] 

+ CF 0 [cos(^ 1 + e t )+ cos(<^ 2 + e 2 ) -2] (15a) 

b 

+ \ l 2( Pi ~ P<Pl + 2 (^2 “ <Pi) 2 ~ P(<P2 " <Pl) 4 ] > 

where 

E = T + + U ( 15b) 

n n n n 

The remainder of the analyses is based on the nondimensional equa- 
tions. In these, the critical buckling load is 

F cr = 1 • < 16 > 

SECTION IV. EQUILIBRIUM 
A. Stable Equilibrium 

At the critical buckling load, equation (4) becomes zero, and the 
necessary conditions for stable equilibrium [ 1] are 
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6*V = 0 


6 4 V> 0 . 


The third variation is identically zero, due to the symmetry of the model. 
Therefore, the fourth variation must be positive definite for stable equilibrium, 
or 


Mr = 


6 a V 


+( 


8*V 

9^1 


T &<Pi + 4 


9 4 V 

d<Pl 9 <P2 


8<pi 6(p 2 + 6 


9*V . , . , 


+ 4 


8 4 V 


dy i 8^2 


7" 5< Pl + 


8 4 V 

9<^2 4 


') 


(17) 


$(p 2 4 > 0 


Using equation ( 15) , setting F = F =1, linearizing at buckling 

cr 

for cpi and <p 2 « 1, and making use of the initial buckling mode of 
equation ( 8) , equation ( 17) is found to be 


24 


<2C - 24(3) 6y? + 480 - 72/3 ( 1 ) 2 dtpf 

•- (18) 

+ 48 ^ (i^I) 3 6c?! 4 + (2C - 1^8 ) ( ^ ^i 4 j > 0 , 


3 ** \f 5* 

where C = . Solution of this shows that the system is in stable 

Li 

equilibrium if 


13 


3 +nT5 
< 12 


or (3 < 0.4363 . . . 


(19) 


The critical value of 0 is therefore defined as 


13 


cr 


3 + n/T 
12 


( 20 ) 
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Because of the symmetry of the model, <5 4 V > 0 as developed in equations 
(18) and (19) meets the sufficiency condition for stable equilibrium. This is 
shown in Reference 4 . 


B. Equilibrium Paths 

To determine the equilibrium paths at bifurcation, equation (14) is 
evaluated at the critical buckling load with an approximation procedure. This 
consists of linearizing in terms of small changes in the variables . Then the 
load change that occurs with small finite displacements of the system can be 
evaluated. Let equation ( 14) be represented by 


f.(<Pi, <p 2 , F) = 0, i = 1,2 . 


With an increment in the variables, the corresponding increment in f. is 


9f. 


af. 


af. 


Af. = 7-1 — <Vi + 6cp 2 + ~r- 6F , 1 = 1,2 . 


i d<Pi 


d(p 2 


8F 


where the derivatives are evaluated in the reference state (<p lf cp 2 , F) 
Since the system remains in equilibrium, 


af, af, 

<Vi + 


af. 


d(pt 


Bn ^ + W 5F = 0 ' i = 1 ' 2 


Substituting equation (14) into this gives 


{2 - 3p[<p 2 + (cp 2 - cpi) 2 ] - CFg cos(<p 1 + e t )} 6<p t 
+ [-1+ 3/3(^ 2 _ ?i) 2 ] <ty 2 ~ [C sin(^ 1 + e^] 6F = 0 , 


(21a) 


[-1+ 3/3(<p 2 “ (Pi ) 2 ] Scpi + [1 “ 3/3(<p 2 " (Pi)* ~ CF cos(<p 2 + e 2 )] 6<p 2 

(21b) 

~[C sin(<p 2 + e 2 )] SF g = 0 . 
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This is placed in the form 


An 8<Pi + A 12 6<P2 = A i3 <5F , 

A 21 6<p i + A 22 5<P2 = A 23 <5F . 

Transposing and using Cramer's rule, this can be solved in terms of either 
6cpi or 8(p 2 ; that is, 

dcp 2 = - ^13 A 22 - A 12 A m ) 6 n > (22a) 

g F = ~_(A u A 22 ^Ai2 A 2l L££l ( ^ 22b ^ 

D = A 13 A 21 - A 11 A 23 . (22c) 


To trace the equilibrium paths, equation (21) was programmed on 
the SDS930 using the form of equation (22) . It is necessary to assign an 
arbitrary or e 2 to make equations (2ia) or (2 lb) nonhomogeneous. 

The program was run for various values of 0 using 8(p t = 0.0001 and 
= 0.01. The results, as shown in Figure 5, indicate that, with /3 < /3 

Cl* 

[equation (20)], the system is in stable equilibrium as the load must increase 

with deflection. With B > B , the load decreases with deflection and the 

cr 

system at buckling is in unstable equilibrium. These correspond to the 

conditions depicted in Figures 3a and 3b, respectively. The unstable system 

(/3 >(3 ) corresponds to the "snap-through" of shells described in Section 

cr 

II. B. 


C. Internal Potential Energy of the System 


The strain energy is a nonlinear function of and ( <p 2 - (p x ) , and 
the total strain energy is the sum of the strain energies of the two springs. 
From equation ( 15) , the nondimensional strain energy is 
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0.02 0.04 0.06 0.08 0.10 


*2 

Figure 5. Equilibrium paths at bifurcation. 

U n = T t2^i 2 + 2 (<^ 2 - <^i) 2 ~jS(<p 2 - «^>i) 4 ] . (23) 

The stationary values for cp^ and ( <p 2 “ <Pt) in equation ( 23) are found by 
taking 

8U 8U 

— — = £ — =0 

d<P 1 9(^2 " Vi) 
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and the maximum value of in each spring is found to occur at 


1 



and cpi - <pi = — — 
N ffi 


( 24 ) 


The energy plot of the strain energy of each spring is shown in Figure 6 . 




Figure 6. Strain energy of the soft spring. 

Since equation ( 24) shows that the maximum value of the strain energy 
is an inverse function of \[]3 , the equilibrium equations are modified by the 
use of 


<p. ' , i 1, 2 , ( 25) 

1 s/jT 

and the equations of equilibrium in the 7). coordinates are 

2r)i - T?i 3 - rj 2 + (rj 2 - r)i) 3 " CF /3 1 / 2 sin = 0 , ( 26a) 

S ^ 
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(26b) 


f )2 - Vi - (V 2 - Vi) 3 - CFg 13'/* sin - 2z - = 0 , 

$ / 


where e. 

i 


has been neglected because 


e i << *i 


The strain energy is 


U n = ~4fT [2V ' " V ' + 2{V2 " Vl) 2 “ (Vz “ V ' } 4] • 


The total spring potential energy of the system is shown in Figure 7 
as a two-dimensional plot of equal energy contours. It can be seen that the 
potential energy is a minimum at the origin where the system is in stable 
equilibrium, and that there are four saddle points at 


(i, 0), (0, 1), (-1, 0), (0, -1), (27) 

which are labeled positions I, II, HI, and IV, respectively, for ready 
reference. The saddle points in the cp coordinates are found by dividing 
these by j3 '/* . 

The potential energy contour plot of Figure 7 is important in the 
analysis of the system. If the system is released from rest at any point 
( or arrives at that point with no kinetie energy) , it will move in the 
direction of a lower potential. If it is outside the saddle point energy contour 
in a region of lower potential, the rotation of the system will increase 
uncontrollably. This is considered to be collapse of the system. 

With a large /3 (very soft spring) , the saddle points are near the 

origin and the system may collapse with small deflections. As /3 decreases, 

the saddle points move out (the actual saddle point for F = 0 is 

b 

(pi = { Pi - cpi) = o.lil . . . radians for j3 = 9. 0, and 1.291. . . radians 
for j3 = 0.6) . At (3 = 0, which corresponds to a linear spring, the saddle 
points are at infinity. 
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Figure 7. Potential energy contour plot of the spring energy. 


D. Change in Saddle Points Under External Loading 

The saddle points determined in Section IV. C were based on the 
internal spring energy only. As an external load is added, the location of 
these saddles will change. At the saddle points for any specified value of 
F , the system can be represented by the equilibrium equations 

s 

1071, %F) = 0, i = i,2 . 
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The incremental change remains in equilibrium, and 


9f. 9f. 9f. 

ST 6 ”* + ** + 77 6F = °' i = 1 - 2 


( 28 ) 


Substitution of equation (26) into equation (28) gives 


2 - 377i 2 - 307 2 - TTi ) 2 - CFg cos 


6vi 


+ [3(t7 2 ~Vi) 2 - 1] <5i7 2 = 


C /S 1 / 2 sin 

e 1 / 2 J 


5 F, 


S ’ 


[3(t7 2 - I7i ) 2 - 1] <5i?i + 


C /3 1 / 2 sin -22- 

^ J 


1 - 3(772 -I7i) 2 -CF,. cos -22- 
s is 1 /* 


6F s ■ 


5i7 2 


where 77 . and F e are evaluated in the reference state. Since these are 
1 b 

linear in 677 ^ 8 r) 2 , and 6 F , they can be solved by Cramer's rule. 

b 

Changing the form to 


Bn <5 t?i + B 12 6772 - B 13 6 Fg 
B21 8 ^ + B 2 2 8 t ]2 = B23 < 5 Fg 


the solution becomes 

6 = (B l3 B 2 2 “B 12 B 2 8 ) 

11 D S’ 
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6t] 2 = i ^ x. g^r_ B - ia.g,2,i) 6Fg ( 

D = B 11 B 22 “ B 12 B2i • 


Substitution of the saddle point values into the above gives the initial 
incremental changes in 7 ) 1 and rj 2 for an incremental external load. The 
initial changes are depicted in Figure 8. 

Figure 8 shows that the changes in the saddle points of positions I 
and III are mirror images about the axes, as are the changes in positions II 
and IV. This is because of the symmetry of the model and the oddness of the 
sine function. The table of Figure 8 also shows that the direction of the saddle 

point movement reverses when p < -jjj, a value less than p^ [equation 
(20)] . Thus, the system becomes unstable under these conditions at a 

n 

different value of /?. The relation between sin — and p is shown in 

n i8 1/2 

Figure 9. 

It is to be expected that the locus of the saddle points, as F c 

s 

increases, will trace a path from at least two of the mirror image positions 
of Figure 8 to the origin. As the saddle point approaches the origin, it is 
also expected that 


lim Fg = F 
Vi, V2^0 

Linearizing equation ( 26) 


= 1 
cr 


for Vi 


and r) 2 near zero gives 


2Vi ~V 2 = CF g Vi , 
V 2 ~ Vi = CF g 7)2 . 
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H «> 9 2 1 .6 ,406 .101 


V 

Figure 9. Relation of sin and j3 . 

( 3 */ ! 


Then as ^ and 0 , 


F , 
S C Cl], • 

(29a) 

F — L.Sl . 

S C Cl], 

(29b) 


Eliminating from equation (29) and solving the resultant quadratic shows 

Vi 

that at r} t = r} 2 = 0 , 


lim 
Vi, t?2 ^ 


F = i = F 
S cr 


E. Locus of Saddle Points 


There are two ways of determining the locus of saddle points as 
F a varies from 0 to 1. One is a direct solution of equation (26) for various 

values of F CT . The other is to eliminate F 0 from equation ( 26) . which 
b b 

gives a single equation represented by 
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ffol, V2, P l/2 ) = 0 


(30) 


where (3 is an arbitrary constant. From this it can be seen that the 
solution is a family of curves for various values of /3. When equation (30) 
has been solved for hi and h 2 > the values can be substituted into equation 
(26) to determine F g , 

The solutions of equation (26) for values of 0 = 0. 6 and (3 = 9.0 
are shown in Figure 10 for the change in saddle point of position I. The change 
in position in is a mirror image of this . The locus of saddle points from 
positions II and IV do not approach the origin because t] 2 is increasing more 
rapidly than r j t . Figure 10 also shows that the solution curves do not vary 
much for widely spread values of j3, but that there is some difference in the 
location of the saddle point for the same values of F Q . 


F. Saddle Point Equal Energy Contours 

As F increases from zero and the saddle points change as shown 
b 

in Figure 10, the entire potential energy contour plot of Figure 7 will change. 
Using equation ( 25) to express the energy in terms of hi and r) 2 yields 


E 

n 


£ 

2/3 


3(hi') 2 + 2rj 1 t r? 2 f cos ^ ' 2 ^ + (r] 2 ') 2 


/ 3 1 / 2 


(31) 


+ CF„ | cos + cos — 


/ 3 1 / 2 


/ 2 2 j +“ [2J]! 2 -hi +2 (t] 2 -rii) 2 


" (h 2 "f?l) 4 ] 


where e has been neglected because £. « h^ . Then the potential energy of 

the external load and the internal springs can be calculated for any load and 

position. Using the values of hi, h 2 » and F that were the basis of Figure 

b 

10, the saddle point equal energy contour lines were calculated for F = 0, 

b 

0. 5, and 0. 9, and are shown in Figure 11 for a j3 = 0. 6. Contour lines for 

saddle points of other values of F would lie in between these and be of the 

b 

same general shape. 
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Figure 10. Saddle point location for values of F 


The change in saddle points of positions II and IV as F increases 

0 25 b 

moves these points outside of the contour V = — — of Figure 7. They 

n jS 1 / 2 

are still saddle points, but at a higher energy level than the saddle points of 
positions I and III. 


The F contour lines of Figure 11 have the same significance as 

D 

those for F 1 = 0 discussed in Section IV. C. The system without kinetic 

s 
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<f>2 -<P 1 



energy will seek a lower potential energy level. If it is in a region outside 
the saddle point energy contour for that F and at a lower potential than the 

s 

saddle point, the system will collapse. 

G. Critical Potential Energy 

Using equation (31) and the values of rj lf 77 2 , and F that determined 

b 

Figure 10, the potential energy of the system can be calculated at the saddle 
points. Then the critical potential energy is 


V 

cr 


= U 


n 



(32) 
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SECTION V. DYNAMICS 


A. Method of Analysis 

The analyses of Section IV were confined to equilibrium states and a 
determination of the critical potential energy for these conditions. This 
section consists of analyses of the dynamic responses to impulsive loads. 
The impulsive loads considered are of two types: 

a. A single square wave such that the impulse, I, is 


I 


F At 
D 


(33) 


where r is the nondimensional time of equation (12). 
b. A pure impulse where 


At = 0 


The pure impulse is an initial velocity of cp x and/or cp 2 from the equilibrium 
state without any load being applied. It represents the limit condition 


I = lim F At 
At-*- 0 D 


The use of a pure impulse is necessary because the results would be influenced 
by the length of At in equation (33). Also the model system of Figure 4 
requires an imperfection, e, to initiate motion when loaded with the impulse 
of equation (33) , while the pure impulse will result in movement of the model 
system without an imperfection. 
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B. Initial Velocities for Impulsive Loads 


When a square -wave impulse is applied to the imperfect system, cp v 
and cp 2 change. The initial velocities of cp± and cp 2 can be determined by 
Lagrange’s impulsive equations [5] . Since the system is at rest prior to the 
impulse, the kinetic energy after the impulse is 


|S— Q. 

9q. i 
1 


(34) 


A 

where Q. is the generalized impulse . Using the nondimens ional energy 
equation (15), equation (34) becomes 


S 2 ^3^’ + $ 2 ' cos( <p 2 + e 2 - <Pi - e t ) J 
= CFAt Sini^i + € t ) 

S 2 [& cos ( (p 2 + e 2 - <P\ - <q) + $2*1 = CFAt Sin (<p 2 + e 2 ) 


where is the initial velocity after the impulse . In the initial condition, 
<Pi = <P 2 = 0* Using equation (33) and linearizing for q and e 2 « 1, the 
initial velocities are found to be 


$i' = ( ei-e 2 ) , (35a) 

Vi' = |pr (3e 2 - €j) , (35b) 


where I = F^Ar 


and S 2 = 


7 - *J~4i 
4 


Therefore, the initial velocities are 


dependent on the imperfections in the system as well as the impulse. This 
also shows, as previously stated, that an imperfection is required for the 
model to move with an impulse as 


$i’ = $ 2 ' = 0 if £i = e 2 = 0 . 
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C. Critical Pure Impulse Velocities 


Section IV. C showed that if the system was at rest outside the outer- 
most saddle point equal energy contour line in a region of lower energy, it 
would collapse. Since the system is conservative, the initial kinetic energy 
is equated to the saddle point potential energy with the system at rest. Then 


E 

n 


= T 


nb 


U .+ 
nf 


a 


nf 


(36) 


where is the total energy, is the initial kinetic energy, and 

and - are the final strain energy and the potential energy of the load, 
respectively, at the saddle point. With equation (32) , equation (36) becomes 


T = U .+ 
nb nf 


0 , = V 
nf cr 


(37) 


and the initial kinetic energy required to reach the saddle point energy is 
equal to the critical potential energy of that saddle point. Since the potential 
of the external load is a negative value [equation ( 15)] , equation (37) shows 
that an external load will decrease the initial kinetic energy required to 
reach the saddle point. 

Equation (37) is evaluated first without an external load. Using 
equation (25) and the energy equation (15) with F = 0, equation (37) is 

O 


s 2 

2/3 


! (v) : 


A A 


+ 2T) 1 '7J 2 , + 




_f_ 
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2 vi - vi+ (tia - ^i ) 2 - (v2 - ’ll ) 4 


(38) 


Transposing and setting and rj 2 at any of the saddle points of Figure 7, 
equation ( 38) becomes 
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2S 


(39) 


3 



+ 


ZriA' +(fi 2 1 ) 2 = 


In the <p coordinate system, this is 


3($ 1 ) 2 + +($ 2 ') 2 



(40) 


Solution of this gives the critical pure impulse velocities for initial 
kinetic energies that will equal the saddle point energy under the conditions 


F 


S 



The solution of equation (39) is the locus of points shown in Figure 13. This 
locus of points is an ellipse with the major axis rotated 22. 5 degrees from 
the vertical. 

Addition of an external load would change the saddle point and reduce 

the energy level on the right side of equation (39) . Thus, the solution for the 

initial pure impulse critical values would be a series of concentric ellipses 

within that shown in Figure 13, with the actual ellipse dependent on F . 

b 

D. Critical Imperfections 


The relationship between the imperfections or eccentricities e* and 

£2 and the critical impulse I to reach the saddle point energy can be 

cr 

determined by substituting equation ( 35) into equation ( 40) to get 


e 2 -2^62+ 3e 2 2 = 


/3C 2 I 2 
cr 


(41) 


where S 2 and C are constants with the values previously determined. The 

solution of equation ( 41) is a locus of points dependent on the value of I . 

cr 
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Figure 14 shows the solution for the value I = 1064. These solutions are 

cr 

also ellipses at an angle of 22. 5 degrees from the coordinate system, and 
have the same eccentricity as that of Figure 13. 

Equation (41) and Figure 14 are based on the condition Fg = 0. If an 
external load is applied to the system, the energy level of the saddle is 
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Figure 14. Critical impulse for imperfections of the system. 

reduced. Thus, the right side of equation (41) would be reduced and the solu- 
tion would be a series of concentric ellipses within that of Figure 14, with the 

actual ellipse dependent on F c . Points outside the ellipse result in increased 

b 

kinetic energy. 

Equation ( 41) gives the values of the imperfections in e. that will 

cause the system to reach the saddle point energy level for specified values 
of the spring softner, ( 3 , and the applied impulse, I. It can be seen that 
as the spring softens {£ increases) or the impulse increases, the critical 
values of the imperfections will decrease, as intuitively expected. To 
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determine the critical impulse to reach the saddle point energy level for 
specified imperfections, equation (41) is rearranged as 


° r Q3 1 / 2 (6 t 2 - 2 €j e 2 + 3c 2 2 ) VZ 

E. Relation of Pure Impulse Velocities to the Imperfect System 

The relation between the pure impulse velocity in Figure 13 and the 
velocities of the imperfect system given by the Lagrange impulsive equations 
can be obtained. For any given conditions of the imperfection and e 2 , 

a relationship is established for <pi' and <p 2 ' by equation (35). This can 
be defined as 


A , „ A . 

<?2 = K . (Pi 


j = 1, 2, 3, 


Converting this by equation ( 25) to 


A A 

Vt' = Krji' 


and substituting into equation ( 39) will give a solution corresponding to a 
point on Figure 13. For example, let e t = 0 and e 2 = £, an arbitrary value. 
Then from equation ( 35) 


A 

<P 1 


t 



4 


and 


A 

<P2 


f 



or 



9 CI t 

3 2S 2 € 



(42) 


Solving equation ( 39) with equation ( 42) will give the values of point A on 
Figure 13. In a similar manner, point B corresponds to e 2 = — e x , and 
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point C to e 2 = ~T € i* Figure 15 shows representative configurations for 

u 

selected points on the ellipse. In this manner, a pure impulse velocity can 
be determined that will represent an imperfect system without the distortion 
introduced in the equations of the imperfect system. 



0 * * 


Figure 15. Configurations represented by pure impulses. 
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Figure 14 is based on the solution of an equation using the imperfec- 
tions and e 2 of the system. It has already been shown that Figure 13 
can be related to imperfections so there is a correspondence between Figures 
13 and 14 and a transformation of coordinates can be established. This 
transformation is 


z i = G l y i 


z 2 “ G 2 J2 


where 


Gi 


G 2 


S 2 

I C /3 1 / 2 
S 2 

I C 0 1 / 2 


2(1- cos 2 ip) 
1 + cos 2 ip 


2(1+ cos 2 ip) 
1 - cos 2 ^ 


i/2 


i/2 


J 


The transformed locations of points A, B, and C of Figure 13 are shown on 
Figure 14. 

F. Responses to Impulsive Loads 


To determine the actual responses to impulsive loads, the nondimen- 

sional equations of motion (13) were programmed for the Univac 1108 using 

a fourth-order Runge-Kutta integration. The different values of F a and /3 

b 

were used with pure impulse velocities determined from equation ( 40) that 
would give the system an initial kinetic energy equal, to the saddle point 
energy. Plots of the responses of the system to these loads and impulses are 
shown in Figures 16 through 21 for a spring softner constant /3 = 0.6. 

The solid lines of Figures 16, 17, and 18 show the response with 
Fg = 0 to the initial velocity of points A, B, and C of Figure 13, respectively. 

In none of these does the configuration reach a point outside the saddle point 
energy contour. 
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0 = 0.6 

Initial kinetic energy to reach saddle potential energy. 

Initial kinetic energy plus 6 percent. 

Figure 16. Response of point A to pure impulse (F = 0) . 

b 
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0 = 0 . 6 

— Initial kinetic energy to reach saddle potential energy. 

— Initial kinetic energy plus 2 percent. 


Figure 17. Response of point B to pure impulse (F = 0) . 
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0 = 0.6 

— Initial kinetic energy to reach saddle potential energy. 
■- Initial kinetic energy plus 16 percent. 


Figure 18. Response of point C to pure impulse ( F 0 = 0) . 

b 
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P = 0.6 

Initial kinetic energy to reach saddle potential energy. 

Initial kinetic energy plus 70 percent. 

Figure 19. Response of point A to pure impulse (F = 0.5). 

b 
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0 = 0 . 6 

Initial kinetic energy to reach saddle potential energy. 

Initial kinetic energy plus 2 percent. 

Figure 20. Response of point B to pure impulse (F c = 0. 5) . 

s 
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<p 2 -<p 1 



j8 =0.6 

— Initial kinetic energy to reach saddLe potential energy. 

— Initial kinetic energy plus 20 percent. 


Figure 2i. Response of point C to pure impulse (F = 0. 5) . 

lb 
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A static load, F = 0.5, was then applied to the system. The 

D 

responses are shown by the solid lines in Figures 19, 20, and 21, which 
correspond respectively to the conditions of Figures 16, 17, and 18, with the 
exception of the added external load. Again, the system does not go outside 
the saddle point energy contour. 

G. Energy Required for Failure 

No equation was found to show the excess energy required to drive 
the system outside the saddle point energy contour and cause the uncontrolled 
deflection and failure of the system previously discussed. Computer runs 
were made with increasing initial kinetic energies until failure of the system 
resulted on the first cycle. The responses are shown by dotted lines in 
Figures 16 through 21. 

For the unloaded configuration, the excess energy required for 

failure varied from less than 2 percent for point B of Figure 14 to almost 

16 percent for point C. For the system loaded at F = 0. 5, the excess 

b 

energy varied from less than 2 percent for point B to nearly 70 percent for 
point A. The excess energy required to make the system fail appears to be 
related to the direction the system initially moves, being lower when it moves 
directly toward a saddle. 

Several of the points of extreme deflection for the systems with excess 
energy in Figures 16 through 21 are outside the saddle point energy contour, 
but in regions of higher energy. Failure only occurs when it is outside this 
contour in a region of lower energy. Because of computer time limitations, 
none of these systems was run long. Without dissipation, it is probable that 
if they ran long enough, all would get in a position where they would fail. 

Figure 16 shows an example of this, as the system did not go over the saddle 
on the approach to + (<p 2 ~ (P\) , but did at -(<p 2 - (px ) . Failure of the system 
in Figure 19 is over the saddle point of position II, which is at a higher 
energy level than position I. A longer run at a lower energy level might put 
it over the saddles shown on the contour plot. However, since all real systems 
have some dissipation energy, early failure is considered a better measure of 
stability. 
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SECTION VI. CONCLUSIONS 


These analyses were conducted with the nondimensional parameters 
P, F, and I to determine their effect on the system, and the major 

conclusions are contained in Paragraphs A through D below. The dimensional 
parameters k, 1, and m were not explicitly considered; however, some 
general comments concerning them are contained in Paragraph E . Areas 
for additional study are listed in Paragraph F. 


A. Effect Of 3 

The spring softening term, p, affects the response in three ways: 


1 . Under static loading, values of p greater than the critical value 

^cr = ’ e< 3 ua ti° n ( 20)J result in unstable equilibrium at the critical 

buckling load and cause a "snap-through" effect as shown in Figure 5. 
However, the unloaded deflected system becomes unstable at the saddle 

point at a different value of pip < ; Fig, 8 

\ 7 r 



2. The location of the saddle point is an inverse function of / 3 1 / 2 

[equation (24)], and increasing p moves the saddle point toward the origin 

(Fig. 7) , which reduces the amount the system can deflect prior to collapse. 

Also, an increase in p decreases the energy of the deflected springs 

[equation (23)] and the critical potential energy (Fig. 12). Since V is 

cr 

equal to the energy required to reach the saddle point, an increase in p 
reduces the dynamic load required to cause collapse. 


3. Changing the value of [3 changes the equation of motion: an 

increase in (3 results in greater amplitudes of the two angles and an increase 

in the period of the system. No analysis of this was made, but the different 

responses for two inputs, identical except for p, are shown in Figure 22. 

p is not important in the dynamic response as there is little difference in 
cr 

responses for values above and below p . 
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B. Effect of the Imperfection € 

1. For the pure impulse (initial velocities of <p^ and <p 2 * equivalent 
to an impulsive load) , the imperfections are not used. However, the initial 
pure impulse velocities correspond directly to particular imperfections 
(Section V. E and Fig. 15) . 

2. The values of ej and e 2 that combine with an impulse load 

to equal the critical potential energy, V , and cause collapse are given by 

the ellipse of Figure 14. For a given impulse, I , all values of and e 2 

on the ellipse result in the same initial kinetic energy. Values of e. within 

the ellipse will give less energy, and values outside the ellipse will give 
greater energy and result in collapse of the system. 

C. Effect of Static Loading 

Increasing the static loading, F_, decreases the kinetic energy that 

D 

the system can absorb without collapse [equation (32) and Fig. 12] . This 
results in smaller concentric ellipses in Figure 13 and reduces the energy 
pulse required to reach the saddle point energy. Loadings greater than 
Fg = 1 will result in collapse without any initial kinetic energy (provided the 
model is perturbed slightly from the origin, <p 1 = cp 2 = 0) . 

D. Effect of Impulses 

An initial pulse of energy in excess of V cr will result in the total 

energy of the system exceeding the saddle point energy, and the system will 
collapse if it reaches a configuration outside the energy contour of the saddle 
point and in a region of lower potential energy ( Fig. 7) . The excess energy 
required for collapse on the first cycle ( a cycle of the system is considered 
to be tne closed path of cp^) varied from less than 2 percent to 70 percent 
and is dependent on the initial conditions. As shown in Figures 16 through 
21, it is less when the system initially moves directly toward the saddle point. 

If a static load is applied, a smaller concentric ellipse would give the velocities; 
however, with a static load, the direction to the saddle changes. 
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Figures 16 through 21 show that the response of the system to an 

impulse is an oscillation of the upper angle ( cp 2 - <p{) in relation to the 

lower angle ((p{), while the latter oscillates more slowly. The rapidity 

of motion of cp 2 - (p^ in relation to cp lt is dependent on the initial 

conditions . Since the response of the system in a configuration space does 

not trace the same path on each cycle, it is probable that, without a 

dissipation function, any system with an initial energy in excess of V will 

cr 

eventually reach a point where it will collapse. 

E. Comments on the Dimensional Parameters 


The effect on the system of parameters k, 1, and m was not 
analyzed. However, certain comments can be made concerning them. These 
parameters are contained in the coefficients of cp 1 and <p 2 in the nondimen- 
sional equations and in the term S 2 [equation (12a)] . 


Changing the assumed mass distribution or making the bars of unequal 
length (Fig. 4) will change the coefficients of cp t ' and cp 2 ' inequations (35), 
(40) , and (41) . This would change the shape and inclination of the response 
ellipses in Figures 13 and 14, resulting in different values for the critical 
impulse. Unequal bar lengths would also change V by changing the relation 

C I* 

between Q, and U. 


In the nondimensional equations, S 2 = 
lowest natural frequency of small oscillations, 
as the lowest natural frequency assumes that ■ 


ml 2 co 2 


where oj is the 

7 - n/U 


The use of S 2 
ml 2 


= 1. Actual values of 


k, 1, and m could change the value of S 2 . Assuming the same relation- 
ship for m l5 m 2 , 1 1# and 1 2 , a change in S 2 would only affect the size 
of the ellipses in Figures 13 and 14. An increase in S 2 would decrease the 
ellipse in Figure 13 and result in smaller pure impulse velocities for collapse. 
This is the same effect as that of an increase in (3. It would increase the size 
of the ellipse in Figure 14 and permit larger imperfections for a given impulse 
without collapse. This is the opposite effect as that of an increase in /3. In 
other words, a system with a higher natural frequency of small oscillations 
(stiffer against small displacements) admits smaller impulsive loading, but 
larger imperfections for a given impulse without collapse. The response for 
three different values of S 2 is shown in Figure 23. 
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The dimensional parameters discussed in this paragraph only affect 
the quantitative values of the analyses and do not change the conclusions. 

F. Areas for Additional Study 

The following are areas where additional analyses would be desirable: 

1. The applicability of the nonlinear equations to very small deflec- 
tions. This includes the effect of 0 on collapse when the saddle point is near 
the origin. 


2. Further analysis of the effect of 0 on the system. 

3. The effects of S 2 and the different dimensional parameters on 
the system. 

4. The effect of the various parameters on 0 . 

5. The applicability of phase -plane trajectories to the analysis. 

6. The excess energy requirements for failure after a number of 
cycles and the effect of the direction of initial displacements. 

7. The extension of this concept to multi-degree -of -freedom systems. 

8. Introduction of dissipation forces. 

G. Applicability to Physical Systems 

These analyses are initial theoretical investigations of a simplified system to 
determine methods that can be used in stability analyses . Further studies as 
outlined above are needed to determine the effectiveness of the methods as the 
system deflections approach linearity or the system becomes more complicated 
(i.e. , more degrees of freedom) . Because of the assumption of frictionless 
pins and the use of nondimensional quantities, direct application of the com- 
pleted analyses to physical systems is not possible. Introduction of dissipation 
forces at the pins and conversion of the equations to dimensional form ( or 
correlation of dimensional characteristics with the numerical values of the 
nondimensional equations as discussed in paragraph VIE) is required to make 
the analyses applicable to physical systems. 
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